(19) 



J 



Europaisches Patentamt 
European Patent Office 
Office europeen des brevets 



III 



(12) 



(43) Veroffentlichungstag: 

09.06.2004 Patentblatt 2004/24 

(21) Anmeldenummer: 03026681.1 

(22) Anmeldetag: 20.11.2003 



(H) EP 1 426 763 A1 

EUROPAISCHE PATENTANMELDUNG 

(51) mt CI7: G01N 33/48, G06F 19/00 



(84) Benannte Vertragsstaaten: 


(71) Anmelder: Bayer Technology Services GmbH 


AT BE BG CH CY CZ DE DK EE ES Fl FR GB GR 


51368 Leverkusen (DE) 


HU IE IT LI LU MC NL PT RO SE SI SK TR 




Benannte Erstreckungsstaaten: 


(72) Erfinder: 


AL LT LV MK 


• Willmann, Stefan, Dr. 




40589 Dusseldorf (DE) 


(30) Prioritat: 03.12.2002 DE 10256315 


• Schmitt, Walter, Dr. 




41470 Neuss (DE) 



(54) Computersystem und Verfahren zur Berechnung eines pharmakokinetischen Verhaltens 
einer chemischen Substanz in Insekten 



CO 
CD 

CO 
CM 



(57) Die Erfindung betrifft ein Computersystem zur 
Berechnung des pharmakokinetischen Verhaltens einer 
chemischen Substanz in Insekten mit: 

einem physiologisch basierten pharmakokineti- 
schen Simulationsmodel (102) eines Insekts zur 
Vorhersage von Konzentrations-Zeit-Verlaufen der 
chemischen Substanz in Kompartimenten des In- 
sekts, wobei das Simulationsmodell zumindest ei- 
nen von der Substanz abhangigen Parameter auf- 
weist, 

einem Vorhersagemodul (1 1 0) zur Vorhersage des 
zumindest einen Parameters basierend auf einer 
physikochemischen Eigenschaft der Substanz. 
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Beschreibung 

[0001] Die Erfindung betrifft ein Computersystem zur Berechnung des pharmakokinetischen Verhaltens einer che- 
mischen Substanz in Insekten auf der Basis eines physiologischen Modells sowie ein entsprechendes Verfahren und 
5 Computerprogrammprodukt. 

[0002] Aus dem Stand der Technik sind zahlreiche physiologisch-basierte Modelle fur Saugetiere bekannt (z. B. 
Charnick et al., J. Pharmacokin. Biopharm. 23, 21 7 (1 995)). Ferner ist aus dem Stand der Technik ein physiologisch- 
basiertes Modell fur Raupen bekannt (Greenwood et al., Pestic. Sci. 30, 97 (1990)). Hierbei wird das Insekt durch 
Kompartimente beschrieben, die jeweils ein einzelnes Organ des Insekts reprasentieren. Die "Verschaltung" der ein- 
10 zelnen Kompartimente ergibt sich aus der bekannten Physiologie von Raupen. Wesentliche Parameter dieses Modells 
sind die Ratenkoeffizienten fur inter-kompartimentellen Massentransport, die die Geschwindigkeit der Verteilung be- 
stimmen, und die Organ-Verteilungskoeffizienten, die das Konzentrationsverhaltnis zwischen dem jeweiligen Organ 
und der Haemolymphe, die der Blutflussigkeit im Saugetierentspricht, im thermodynamischen Gleichgewichtangeben. 
[0003] Aus dem Stand der Technik ist es bekannt, solche kompartimentellen Modelle einzusetzen, um ein expen- 
ds mentell ermitteltes pharmakokinetisches Profil einer Substanz retrospektiv durch Anpassung der Raten- und Vertei- 
lungskoeffizienten zu beschreiben (z. B. Lagadic et al., Pestic. Biochem. Physiol. 45, 105 (1993) & Pestic. Biochem. 
Physiol. 48, 173 (1994)). 

[0004] Der Erfindung liegt demgegenuber die Aufgabe zugrunde ein Computersystem zur Vorhersage des pharma- 
kokinetischen Verhaltens einer chemischen Substanz in einem Insekt zu schaffen, sowie ein entsprechendes Verfahren 
20 und Computerprogrammprodukt. 

[0005] Die der Erfindung zugrunde liegenden Aufgaben werden mit den Merkmalen der unabhangigen Patentan- 
spruche jeweils gelost. Bevorzugte Ausfuhrungsformen der Erfindung sind in den abhangigen Patentanspruchen an- 
gegeben. 

[0006] Die Erfindung ermoglicht es, auf besonders effiziente Art und Weise eine Vorhersage des pharmakokineti- 
25 schen Verhaltens einer chemischen Substanz in einem Insekt zu berechnen. Insbesondere ermoglicht es die Erfindung, 
die Aufnahme, Verteilung und Ausscheidung von chemischen Substanzen in Insekten auf der Basis physikochemischer 
Parameter abzuschatzen. 

[0007] Hierzu wird ein physiologisch basiertes pharmakokinetisches Simulationsmodell eines Insekts zur Vorhersage 
von Konzentrations-Zeit-Verlaufen einer chemischen Substanz in den Kompartimenten des Insekts verwendet. Das 
30 Simulationsmodell beinhaltet zumindest einen von der zu untersuchenden Substanz abhangigen Parameter. Der oder 
die Parameter des Simulationsmodells werden fur eine bestimmte Substanz basierend auf einer oder mehrerer der 
physikochemischen Eigenschaften der Substanz vorhergesagt. 

[0008] Bei diesen physikochemischen Parametern handelt es sich beispielsweise um die Lipophilie der Substanz 
beschrieben durch den Verteilungskoeffizienten zwischen Wasser und Phospholipidmembranen oder den Oktanol- 
35 Wasser-Verteilungskoeffizienten, das Molekulgewicht oder die Loslichkeit. Die relevanten physikochemischen Para- 
meter der Substanz konnen entweder aus einfachen Experimenten bestimmt oder mittels an sich bekannter Verfahren 
wie QSAR (Quantitative Structure Activity Relations) oder neuronalen Netzen direkt aus dem Deskriptor der chemi- 
schen Struktur der Substanz ermittelt werden. 

[0009] Im letzteren Fall ist es mit dem erfindungsgemaBen Verfahren sogar moglich, virtuelle, d.h. noch nicht syn- 
40 thetisierte Substanzen im Hinblick auf ihre Aufnahme- und Verteilungseigenschaften in Insekten zu bewerten. Weiterhin 
konnen aufgrund des etablierten Zusammenhangs zwischen physikochemischen und pharmakokinetischen Eigen- 
schaften allgemeine Kriterien zur Optimierung von Insektizidenwirkstoffen abgeleitet werden. 

[0010] Nach einer bevorzugten Ausfuhrungsform der Erfindung wird als substanzabhangiger Parameter des Simu- 
lationsmodells der Ratenkoeffizient des inter-kompartimentellen Massentransports verwendet, welcher proportional 

45 zum Produkt aus der Permeabilitatfur die Substanz und der effektive Oberflache der Kompartimente ist. Vorzugsweise 
gibt es also einen substanzabhangigen Parameter fur jedes Kompartiment des Simulationsmodells, der die Permea- 
bilitat des betreffenden Kompartiments fur die Substanz und die effektive Oberflache des betreffenden Kompartiments 
beinhaltet. Der Permeabilitatskoeffizient beschreibt den Substanzfluss durch die Zellmembran. 
[0011] Besonders vorteilhaft ist, dass zur Ermittlung der Parameter fur das Simulationsmodell basierend auf einer 

50 physikochemischen Eigenschaft der Substanz keine weiteren experimentellen Untersuchungen notwendig sind, son- 
dern das der oder die Parameter basierend auf der physikochemischen Eigenschaft der Substanz ermittelt werden 
kann. Diese Ermittlung erfolgt auf der Grundlage einer Datenbasis, die fur verschiedene Tests ubstanzen zuvor expe- 
rimentell ermittelt worden ist. Die Datenbasis beinhaltet die fur die Testsubstanzen experimentell ermittelten substanz- 
abhangigen Parameter des Simulationsmodells sowie die physikochemischen Eigenschaften der Testsubstanzen. Die- 

55 se Datenbasis wird zur Pradiktion des oder der substanzabhangigen Parameter fur eine zu untersuchende Substanz 
verwendet. 

[0012] Nach einer bevorzugten Ausfuhrungsform der Erfindung wird aus der Datenbasis durch eine lineare Regres- 
sion eine Berechnungsfunktion gewonnen. Beispielsweise handelt es sich bei der Berechnungsfunktion um eine Funk- 
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tion der Lipophilie und des Molekulargewichts. Zur Vorhersage eines Parameters fur das Simulationsmodell fur eine 
zu untersuchende Substanz mussen also lediglich die Lipophilie und das Molekulargewicht der zu untersuchenden 
Substanz mit der Berechnungsfunktion ausgewertet werden, um den Parameter zu erhalten. Mit dem vorhergesagten 
Parameter kann dann eine konkrete Simulation von Konzentrations-Zeit-Verlaufen bei der Aufnahme und Ausschei- 
5 dung der Substanz in dem Insekt durchgefuhrt werden. 

[0013] Alternativ konnen anstelle der durch lineare Regression gewonnenen Berechnungsfunktion auch andere an 
sich bekannte Pradiktionsverfahren verwendet werden. 

[0014] Von besonderem Vorteil ist, dass nachdem einmal eine Datenbasis fur eine bestimmte Anzahl von Testsub- 
stanzen ermittelt worden ist, fur die Simulation des pharmakokinetischen Verhaltens weiterer Substanzen keine zu- 
10 satzlichen Experimente mehr erforderlich sind. Dies ermoglicht es mit einem hohen Durchsatz "Kandidaten" fur po- 
tenzielle Insektizide hinsichtlich deren pharmakokinetischen Verhalten zu beurteilen. Dadurch ist eine erhebliche Be- 
schleunigung bei der Auffindung, Entwicklung und Optimierung von neuen Insektiziden moglich. 
[0015] Im Weiteren werden bevorzugte Ausfuhrungsformen der Erfindung mit Bezugnahme auf die Zeichnungen 
naher erlautert. Es zeigen: 

15 

Figur 1 ein Blockdiagramm einer Ausfuhrungsform eines erfindungsgemaBen Computersystems, 



Figur 2 ein Flussdiagramm einer Ausfuhrungsform des erfindungsgemaBen Verfahrens, 

20 Figur 3 ein physiologisch basiertes Modell einer Raupe als Grundlagefur ein erfindungsgemaBes Simulationsmodell 
mit einem substanzabhangigen Parameter, 

Figur 4 ein Diagramm zur Darstellung derfurdas Simulationsmodell der Figur 3 gewonnenen Datenbasis zur Durch- 
fuhrung einer linearen Regression, um eine Berechnungsfunktion zu ermitteln. 

25 

[0016] Die Figur 1 zeigt ein Computersystem 100. Dabei kann es sich um einen ublichen Personalcomputer (PC), 
eine Workstation oder um ein Client-Serversystem handeln. 

[0017] Das Computersystem 100 hat ein Simulationsmodell 102. Hierbei handelt es sich um ein physiologisch ba- 
siertes pharmakokinetisches Modell eines Insekts. Hierzu bildet das Simulationsmodell 102 die Kompartimente des 
30 Insekts ab, und ermoglicht so die Vorhersage von Konzentrations-Zeit-Verlaufen einer Substanz in den Kompartimen- 
ten des Insekts. 

[001 8] Zu dem Simulationsmodell 1 02 gehoren physiologische Parameter, 1 04 die nur von der Art des zu beschrei- 
benden Insekts abhangen, sowie ein oder mehrere substanzabhangige Parameter 106. Zur Durchfuhrung einer Si- 
mulationsberechnung von Konzentrations-Zeit-Verlaufen ist also die Eingabe des Parameterwerts des substanzab- 

35 hangigen Parameters 106 fur die auszuwertende Substanz erforderlich. Ein konkretes Beispiel fur eine Ausfuhrungs- 
form des Simulationsmodells 102 wird weiter unten mit Bezugnahme auf die Figuren 3 und 4 naher erlautert. 
[0019] Das Computersystem 100 hat ferner eine Datenbank 108. Die Datenbank 108 dient zur Speicherung einer 
Datenbasis, die auf der Grundlage von experimentellen Untersuchungen des pharmakokinetischen Verhaltens von 
Tests ubstanzen in dem Insekt gewonnen worden sind. Fur jede derzuvor experimentell untersuchten Testsubstanzen 

40 beinhaltet die Datenbank 1 08 den oder die Parameterwerte, die fur die betreffende Testsubstanz experimentell ermittelt 
worden sind sowie zumindest eine physikochemische Eigenschaft der betreffenden Testsubstanz. 
[0020] Diese in der Datenbank 108 gespeicherte Datenbasis bildet die Grundlage fur die Vorhersage eines Para- 
meterwerts fur eine neue zu untersuchende Substanz in dem Pradiktionsmodul 110. Beispielsweise wird durch ein 
lineares Regressionsverfahren aus der in der Datenbank 108 gespeicherten Datenbasis eine Berechnungsvorschrift 

45 gewonnen, die es ermoglicht, aus einer physikochemischen Eigenschaft einer zu untersuchenden Substanz, wie bei- 
spielsweise deren Lipophilie oder deren Molekulargewicht, den fur die Durchfuhrung der Simulation mit dem Simula- 
tionsmodell 102 erforderlichen substanzabhangigen Parameterwert zu gewinnen. 

[0021 ] Das Computersystem 1 00 hat ferner ein Eingabe-/Ausgabe-Modul 1 1 2 zur Eingabe der physikochemischen 
Eigenschaft einer zu untersuchenden Substanz. Ferner erfolgt uber das Eingabe-/Ausgabe-Modul 112 die Ausgabe 

50 der simulierten Konzentrations-Zeit-Verlaufe. 

[0022] Das Eingabe-/Ausgabe-Modul 112 kann mit einer Datenbank gekoppelt sein, in der sich Deskriptoren von 
real existenten oder virtuellen Substanzen, wie z.B. potenziellen Insektiziden, und deren physikochemische Eigen- 
schaften befinden. Bei noch nicht synthetisierten, virtuellen Substanzen kann die fur die Eingabe in das Computersy- 
stem 100 erforderliche physikochemische Eigenschaft direkt aus dem Deskriptor der chemischen Struktur der Test- 

55 substanz mittelsan sich bekannter Verfahren wie QSAR oder neuronalen Netzen ermittelt werden. Im Fall von virtuellen 
Substanzen kann anstelle der physikochemischen Eigenschaft auch der Deskriptor uber das EingabeVAusgabe-Modul 
1 1 2 eingegeben werden. Auf der Grundlage des Deskriptors wird dann in dem Computersystem 1 00 selbst die physi- 
kochemische Eigenschaft als EingangsgroBe fur das Pradiktions-Modul 110 ermittelt. 
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[0023] Die Simulationsergebnissewerden uber das Eingabe-/Ausgabe-Modul 1 12z.B. in die Datenbankeingegeben, 
sodass zu einem spateren Zeitpunkt eine Bewertung der Simulationsergebnisse erfolgen kann. 
[0024] Die Figur 2 zeigt ein entsprechendes Flussdiagramm. In dem Schritt 200 wird eine physikochemische Eigen- 
schaft einer zu untersuchenden Substanz eingegeben. Bei der physikochemischen Eigenschaft handelt es sich bei- 

5 spielsweise um die Lipophilie.oder das Molekulargewicht der Substanz. 

[0025] In dem Schritt 202 werden die Parameterwerte der substanzabhangigen Parameter des Simulationsmodells 
fur die Substanz auf der Grundlage einer Datenbasis basierend auf der physikochemische Eigenschaft der zu beur- 
teilenden Substanz berechnet. Diese Datenbasis beinhaltet zuvor experimentell ermittelte Parameterwerte von ver- 
schiedenen Testsubstanzen. Diese Berechnung erfolgt in dem Schritt 202. 

10 [0026] In dem Schritt 204 werden die in dem Schritt 202 berechneten Parameter in das Simulationsmodell eingege- 
ben. Dort erfolgt dann die Berechnung von Konzentrations-Zeit-Verlaufen der zu untersuchenden Substanz in den 
Kompartimenten des Insekts. In dem Schritt 206 werden diese Konzentrations-Zeit-Verlaufe ausgegeben und konnen 
dann bewertet werden. 

[0027] Zur Ermittlung einer moglichst aussagekraftigen Datenbasis fur die Pradiktion der Parameterwerte ist es vor- 
15 teilhaft, wenn die zur experimentellen Ermittlung der Datenbasis verwendeten Testsubstanzen chemisch moglichst 
divers sind. Vorzugsweise sollte eine neue zu untersuchende Substanz, fur die die Parameterwerte vorhergesagt wer- 
den sollen, innerhalb des von den Testsubstanzen reprasentierten Versuchsraums liegen. 

[0028] Die Figur 3 zeigt beispielhaft ein physiologisches Modell fur eine Raupe. Dieses physiologische Modell basiert 
auf dem von Greenwood et al. (siehe oben) angegebenen Modell. 
20 [0029] Zu dem Modell gehoren die folgenden acht Kompartimente der Raupe: Haemolymphe (hi) 300, Fett-Korper 
(fb) 302, Muskel (mu) 304, Oberflache der Kutikula (cs) 306, Kutikula (c) 308, Darmwand (gw) 310, Darminhalt (gc) 
312 und Nervenstrang (mc) 314. Eine zu untersuchende Substanz wird entweder topisch uber die Kutikula 308, oral 
uber den Darm oder durch direkte Injektion in die Haemolymphe 300 verabreicht. 

[0030] Die Haemolymphe 300 dient, ahnlich wie das Blut bei Saugetieren, als Haupttransportphase zwischen den 
25 verschiedenen Organen. 

[0031 ] Stoffwechsel f indet mit der Rate kc in der Kutikula 308, mit der Rate khl in der Haemolymphe 300 und mit der 
Rate kgw in der Darmwand 31 0 statt. Die Ausscheidung f indet mit der Rate kg C statt. 

[0032] Die Haemolymphe 300 wird als frei zirkulierende Fliissigkeit modelliert, welche in Kontakt mit den durch die 
Kompartimente beschriebenen Raupenorganen steht. Transportvorgange zwischen den Kompartimenten erfolgen 
30 durch passive Diffusion mit Permeation uber die Membranen als ratenlimitierenden Schritt. Die Ratenkoeffizienten \ 
des Massentransports zwischen den Kompartimenten werden durch das Permeabilitats-Oberflachenprodukt P X A X und 
das Volumen V x der Organe bestimmt: 

P A 

35 K = -rr^ ( x G {c,mu,fb,nc,gw}) (1) 

V x 

[0033] Hierin bedeutet 
40 p x die Permeabilitat der Membran eines Organes x fur die Substanz, 
A x die Oberflache der Menbran des Organes x 
c Kutikula 

45 

mn Muskel 
fb Fettkorper 
50 nc Nervenstrang 
gw Darmwand 

[0034] Der Gleichgewichtszustand wird dabei nach einer Zeit t x » Vy (P X A X ). erreicht. Das Verhaltnis der Konzen- 
55 trationen in den peripheren Kompartimenten im Gleichgewichtszustand wird dabei durch die Verteilungskoeffizienten 
bezogen auf die Haemolymphe K x := K x/N : bestimmt: 
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K x = pJL (x G {c,mu,fb,nc,gw}) (2) 

5 [0035] Bei K x handelt es sich also um den Verteilungskoeffizienten einer Substanz zwischen Haemolymphe und dem 
Organ x im Gleichgewichtszustand. 

[0036] Zusatzlich mussen die Verteilungskoeffizienten zwischen der Oberflache der Kutikula 306 und der Kutikula 
(K c/cs ) sowie auch der Verteilungskoeffizient zwischen der Darmwand 310 und dem Darminhalt 312 (K gc/gw ) bekannt 
sein. 

10 [0037] Das von Greenwood et al. bekannte Modell mit acht Kompartimenten ist zur Berucksichtigung von 
verschiedenen Formen der Verabreichung (oral, topisch oder durch Injektion in die Haemolymphe), 
einen 11 Open Loop"-Massentransport von Nahrung Liber den Ernahrungskanal und 

15 

die Zeitabhangigkeit der Organvolumina wegen des Larvenwachstums modifiziert worden. 

[0038] Basierend auf diesem biophysikalischen Modell der Figur 3 lasst sich fur jedes Organ x der Raupe eine Mas- 
sengleichgewichtsbeziehung durch eine Differenzialgleichung beschreiben. Solche Massengleichgewichtsbeziehun- 
20 gen sind fur die Organe der Raupe im Anhang angegeben (Formeln A1 bis A8). 

[0039] Zur Berechnung eines pharmakokinetischen Profils mittels der Massengleichgewichtsbeziehungen der Or- 
gane der Raupe mussen insgesamt 23 Parameterwerte bekannt sein. Hierbei handelt es sich um 

die Organvolumina V x (acht Parameterwerte sowie zusatzliche Parameterwerte zur Beschreibung deren zeitlicher 
25 Veranderungen aufgrund des Larvenwachstums), 

die Verteilungskoeffizienten K x der Subtanz zwischen Haemolymphe 300 und peripheren Kompartimenten (c, mu , 
fb, nc und gw), 

30 - die Verteilungskoeffizienten zwischen der Oberflache der Kutikula und der Kutikula sowie zwischen der Darmwand 
und dem Darminhalt, 

die Stoffwechselratenkonstanten in der Kutikula, der Haemolymphe und der Darmwand, sofern die zu untersu- 
chende Substanz metabolisiert wird, und 

35 

Ratenkoeffizienten X x = [P X A]^\/ X fur den Massentransport zwischen den Kompartimenten (c, mu, fb, nc und gw). 

[0040] Die Verteilungskoeffizienten sowie die Permeabilitats-Oberflachenprodukte und damit die Ratenkoeffizienten 
hangen von physiologischen Parametern des Insekts ab sowie von den physikochemischen Eigenschaften der zu 
40 untersuchenden Substanz. Fernersind auch die Stoffwechselratenkonstanten substanzspezifisch. 

[0041] Die Organvolumina konnen durch an sich bekannte Methoden experimentell ermittelt werden. Die Vertei- 
lungskoeffizienten K x der Substanz zwischen Haemolymphe und Organ x konnen wirfolgt berechnet werden: 



X ~ K h | /H o Kx / H 2° ~ f water,x + K fat f fat,x + K protein f prvtein,x ( 3 ) 



Verteilungskoeffizient der Substanz im Gleichgewicht zwischen Wasser und Fett (Lipophilie) 
Verteilungskoeffizient der Substanz im Gleichgewicht zwischen Wasser und Protein 
Volumenanteil von Wasser in dem Kompartiment x, 
Volumenanteil von Fett in dem Kompartiment x, 
Volumenanteil von Protein in dem Kompartiment x. 



wobei 
50 K tat = 

^protein = 
A/vater.x = 

55 

f fat,x = 
^protein, x = 
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[0042] Die Membranaffinitat (MA) oder alternativ der Oktanol/Wasser-Verteilungskoeffizient (K o/w ) der Substanz 
kann als Schatzung des Verteilungskoeffizienten K fat verwendet werden. Der Verteilungskoeffizient K protejn kann bei- 
spielsweise aus der Bindungskonstanten von Humanserumalbumin (K d HSA in [mmol]) und dem Molekulargewicht des 
Proteins (65 kDa) ermittelt werden: 

5 

hsa =;t4a /65 W 
K d 

10 [0043] Die Formel 3 kann durch ^Combination der Fett- und Proteinanteile in einem einzigen organischen Anteil ver- 
einfacht werden: 

K x/h 2 o 

K x = i< mit K x/H 9 0 = f water,x + 0 " f water.x) K fat ( 5 ) 

15 ^hl/H 2 0 2 

[0044] Die Ratenkonstanten des Massentransports zwischen Kompartimenten werden experimentell fur die Test- 
substanzen ermittelt, indem sie mit experimentellen pharmakokinetischen Daten abgeglichen und mit den physikoche- 
mischen Eigenschaften der Testsubstanzen korreliert werden. 
20 [0045] Approximativ kann angenommen werden, dass der Permeabilitatskoeffizient P fur die Substanz fur alle Or- 
gane x identisch ist, d. h. Px = P fur alle Organe x, mit Ausnahme der Kutikula. Ferner ist P proportional zur Lipophilie 
(Kfat= MA oder K fat = K o/w , d. h. dem Verteilungskoeffizienten Oktanol/Wasser) und dem Membrandiffusionskoeffizi- 
enten D mem der Substanz: 

P - K fat D mem (6) 

[0046] Mittels einer Exponentialbeziehung wird die Abhangigkeit des Membrandiffusionskoeffizienten von dem Mo- 
lekulargewicht (MW) der Substanz beschrieben: 

30 

D mem - MW" Smem (7) 
[0047] Die Formeln 6 und 7 konnen miteinander kombiniert werden, so dass 

35 

P = a K fat MW _Smem (8) 

resultiert, wobei a eine Konstante ist. Mittels Gleichung 1 konnen die Ratenkonstanten des Massentransports zwischen 
den Kompartimenten wie folgt ausgedruckt werden: 

^x = « J* ^ MW _Smem (9) 
v x 

45 

[0048] Durch Logarithmierung erhalt man daraus 

Log(^) = Log(a^)-s mem Log(MW) (1 0) 

50 Nat V x 

[0049] Die Figur 4 zeigt eine entsprechende doppelt logarithmische Darstellung von X^K fat uber dem Molekularge- 
wicht MW mehrerer Testsubstanzen fur verschiedene Kompartimente. Ebenso zeigt die Figur 4 durch lineare Regres- 
sionen ermittelte Ausgleichsgeraden durch die Messpunkte. 

[0050] Mittels linearer Regression lasst sich also aus der experimentell ermittelten Datenbasis die Steigung s mem 
bestimmen sowie auch der Achsenabschnitt otA x /V x . Damit hat man eine Berechnungsvorschrift erhalten, die es er- 
laubt, allein aufgrund der Lipophilie (K fat ) und des Molekulargewichts MW einer neuen, auch virtuellen Substanz die 
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GroBe X x zu berechnen. Aufgrund der Formeln 8 und 9 erhalt man so P • A x . Auf der Basis dieses substanzabhangigen 
Parameters sowie der Gleichgewichts-Verteilungskoeffizienten K x bzw. und K gc/gw aus Gl. (3) bzw. (5) lasst sich 
so das Gleichungssystem im Anhang z. B. numerisch losen, und man erhalt die Zeit-Konzentrations-Verlaufe der zu 
untersuchenden Substanz in den einzelnen Kompartimenten. 

5 

Anhang 

[0051] Gleichungssystem fur das pharmakokinetische Simulationsmodell der Raupe 
10 (cs) 



15 



25 



30 



35 



40 



45 



50 



V cs ^ = -[PA] c (C cs - jr ^) (A1) 



(C) 



V c ^ = [PA] C (C cs - + [PA] C (C h1 - % - k c C c (A2) 



(M) 



dC c 
Vhl = ~ Zt PA ]x (C hl - -2-) - k hl C h , (x € {c,mu,fb,nc,gw}) (A3) 



(mu) 



(fb) 



V mu ^p = [PA] mu (C h1 -^) (A4) 



V fb ^ = [PA] fb (C h1 - R ^) (A5) 



(nc) 



dC„ 



V nc ^ = [PA] nc (C h1 -^ c ) (A6) 



(gw) 



WL [X nc 



V gw - Hp - [ PA lgw ( C h1 " j<^) " [ PA ]gw ( C gw ' « 9 ° ) ' k gw C gw ( A7 ) 
UL ,x gw rx gc/gw 
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v gc^ = r pA W c gw-j^) (A8) 

rx gc/gw 

5 C x : Konzentration der Substanz in dem Kompartiment x, 

V x : Volumen des Kompartiments x; das Volumen kann nach einer experimentell ermittelten Funktion zeitlich 

veranderlich sein (V x (t)), 

[PA] X : Permeabilitat-Oberflachenprodukt von Kompartiment x, 

K x/y : Verteilungskoeffizient zwischen den Kompartimenten x und y, 

10 K x = K x/h) : Verteilungskoeffizient zwischen Organ x und Haemolymphe, 

k x : Ratenkonstante fur den Stoffwechsel (x = c, hi, gw) 
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35 

1 . Computersystem zur Berechnung eines pharmakokinetischen Verhaltens einer chemischen Substanz in Insekten 
mit: 

einem physiologisch basierten pharmakokinetischen Simulationsmodel (102) eines Insekts zur Vorhersage 
40 von Konzentrations-Zeit-Verlaufen der chemischen Substanz in Kompartimenten des Insekts, wobei das Si- 

mulationsmodell zumindest einen von der Substanz abhangigen Parameter aufweist, 

einem Vorhersagemodul (11 0) zur Vorhersage des zumindest einen Parameters basierend auf einer physiko- 
chemischen Eigenschaft der Substanz. 

45 

2. Computersystem nach Anspruch 1 , wobei es sich bei dem Parameter urn das Produkt aus der Permeabilitat fur 
die Substanz und der effektiven Oberflache der Kompartimente handelt. 

3. Computersystem nach Anspruch 1 oder 2, wobei es sich bei dem Parameter urn den Ratenkoeffizienten fur inter- 
so kompartimentellen Massentransport (X x = P X A X /V X ) handelt, wobei insbesondere das Volumen zumindest eines 

der Organe des Insekts eine Funktion derZeit ist (V x =V x (t)). 

4. Computersystem nach Anspruch 1 , 2 oder 3, wobei es sich bei dem Parameter urn den Gleichgewichts-Vertei- 
lungskoeffizienten zwischen Organ und Haemolymphe (K x , x e {c, mu, fb, nc und gw}), zwischen Oberflache der 

55 Kutikula und Kutikula (K^), oder zwischen Darmwand und Darminhalt (Kg C/gw ) handelt. 

5. Computersystem nach einem der vorhergehenden Anspruche 1 bis 4, wobei es sich bei der physikochemischen 
Eigenschaft urn den Verteilungskoeffizienten zwischen Wasser und Phospholipidmembranen, den Oktanol-Was- 



8 



EP 1 426 763 A1 



ser-Verteilungskoeffizienten, das Molekulgewicht, die Loslichkeit, und/oder eine Kombination dieser Parameter 
der Substanz handelt. 

6. Computersystem nach einem der vorhergehenden Anspruche 1 bis 5, mit einem QSAR-Modul oder einem neu- 
5 ronalen Netz zur Bestimmung der physikochemischen Eigenschaft aus einem Deskriptor der chemischen Struktur 

der Substanz. 

7. Computersystem nach einem der vorhergehenden Anspruche 1 bis 6, wobei das Vorhersagemodul auf einer Da- 
tenbasis (108) basiert, die die physikochemischen Eigenschaften von Testsubstanzen und experimentell fur die 

10 Testsubstanzen ermittelten Parameter beinhaltet. 

8. Computersystem nach einem der vorhergehenden Anspruche 1 bis 7, wobei das Vorhersagemodul eine Berech- 
nungsfunktion zur Berechnung des Parameters aus der Lipophilic und/oder dem Molekulargewicht der Substanz 
beinhaltet. 

15 

9. Computersystem nach Anspruch 8, wobei die Berechnungsfunktion auf einer linearen Regression experimentell 
ermittelter Parameterwerte basiert. 

10. Verfahren zur Berechnung eines pharmakokinetischen Verhaltens einer chemischen Substanz mit einem physio- 
20 logisch basierten pharmakokinetischen Simulationsmodell eines Insekts zur Vorhersage von Konzentrations-Zeit- 

Verlaufen der chemischen Substanz in Kompartimenten des Insekts, wobei das Simulationsmodell zumindestens 
einen von der Substanz abhangigen Parameter aufweist, mitfolgenden Schritten: 

Eingabe einer physikochemischen Eigenschaft der Substanz in ein Vorhersagemodul zur Vorhersage des 
25 zumindest einen Parameters fur die Substanz, 

Durchfuhrung einer Simulation mit Hilfe des Simulationsmodells zur Vorhersage von Konzentrations-Zeit-Ver- 
laufen der chemischen Substanz basierend auf dem vorhergesagten zumindest einen Parameter. 

30 11. Verfahren nach Anspruch 10, wobei es sich bei dem Parameter urn das Produkt aus der Permeabilitat fur die 
Substanz und der effektiven Oberflache der Kompartimente handelt. 

12. Verfahren nach Anspruch 10 oder 11, wobei es sich bei dem Parameter urn den Ratenkoeffizienten fur inter- 
kompartimentellen Massentransport ( x = P X A X /V X ) handelt, wobei insbesondere das Volumen zumindest eines der 

35 Organe des Insekts eine Funktion der Zeit ist (V x = V x (t)). 

13. Verfahren nach Anspruch 10, 11 oder 12, wobei es sich bei dem Parameter urn den Gleichgewichts-Verteilungs- 
koeffizienten zwischen Organ und Haemolymphe (K x , x £ {c, mu, fb, nc und gw}), zwischen Oberflache der Kutikula 
und Kutikula (K c/cs ), oder zwischen Darmwand und Darminhalt (Kg C/gw ) handelt. 

40 

14. Verfahren nach einem der vorhergehenden Anspruche 10 bis 13, wobei es sich bei der physikochemischen Ei- 
genschaft urn den Verteilungskoeffizienten zwischen Wasser und Phospholipidmembranen, den Oktanol-Wasser- 
Verteilungskoeffizienten, das Molekulgewicht, die Loslichkeit, und/oder eine Kombination dieser Parameter der 
Substanz handelt. 

45 

15. Verfahren nach einem der vorhergehenden Anspruche 10 bis 14, wobei die physikochemische Eigenschaft der 
Substanz mittels eines QSAR oder mittels einem neuronalen Netz bestimmt wird. 

16. Verfahren nach einem der vorhergehenden Anspruche 10 bis 15, wobei die Vorhersage des zumindest einen 
50 Parameters basierend auf physikochemischen Eigenschaften von Testsubstanzen und experimentell fur die Test- 
substanzen ermittelten Parameterwerten erfolgt. 

17. Verfahren nach einem der vorhergehenden Anspruche 10 bis 16, wobei die Vorhersage des zumindest einen 
Parameters mit einer Berechnungsfunktion aus der Lipophilie und/oder dem Molekulargewicht der Substanz er- 

55 folgt. 

18. Verfahren nach einem der vorhergehenden Anspruche 1 0 bis 1 7, wobei die Berechnungsfunktion auf einer linearen 
Regression experimentell ermittelter Parameterwerte basiert. 
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19. Computerprogrammprodukt, insbesondere digitales Speichermedium, mit Programmmitteln zur Berechnung eines 
pharmakokinetischen Verhaltens einer chemischen Substanz mit einem physiologisch basierten pharmakokineti- 
schen Simulationsmodell eines Insekts zur Vorhersage von Konzentrations-Zeit-Verlaufen der chemischen Sub- 
stanz in Kompartimenten des Insekts, wobei das Simulationsmodell zumindestens einen von der Substanz ab- 

5 hangigen Parameter aufweist, mitfolgenden Schritten: 

Eingabe einer physikochemischen Eigenschaft der Substanz in ein Vorhersagemodul zur Vorhersage des 
zumindest einen Parameters fur die Substanz, 

10 - Durchfuhrung einer Simulation mit Hilfe des Simulationsmodells zur Vorhersage von Konzentrations-Zeit-Ver- 

laufen der chemischen Substanz basierend auf dem vorhergesagten zumindest einen Parameter. 

20. Computerprogrammprodukt nach Anspruch 15, wobei es sich bei dem Parameter um das Produkt aus der Per- 
meabilitat fur die Substanz und der effektiven Oberflache der Kompartimente handelt. 

15 

21. Computerprogrammprodukt nach Anspruch 19 oder 20, wobei es sich bei dem Parameter um den Ratenkoeffizi- 
enten fur inter-kompartimentellen Massentransport (X x = P X A X /V X ) handelt, wobei insbesondere das Volumen zu- 
mindest eines der Organe des Insekts eine Funktion der Zeit ist (V x = V x (t)). 

20 22. Computerprogrammprodukt nach Anspruch 19, 20 oder 21 , wobei es sich bei dem Parameter um den Gleichge- 
wichts-Verteilungskoeffizienten zwischen Organ und Haemolymphe (K x , x e {c, mu, fb, nc und gw}), zwischen 
Oberflache der Kutikula und Kutikula (K c/cs ), oder zwischen Darmwand und Darminhalt (K gc/gw ) handelt. 

23. Computerprogrammprodukt nach einem der vorhergehenden Anspruche 19 bis 22, wobei es sich bei der physi- 
25 kochemischen Eigenschaft um den Verteilungskoeffizienten zwischen Wasser und Phospholipidmembranen, den 

Oktanol-Wasser-Verteilungskoeffizienten, das Molekulgewicht, die Loslichkeit, und/oder eine Kombination dieser 
Parameter der Substanz handelt. 

24. Computerprogrammprodukt nach einem der vorhergehenden Anspruche 1 9 bis 23, mit einem QSAR-Modul, oder 
30 einem neuronalen Netz zur Bestimmung der physikochemischen Eigenschaft aus einem Deskriptor der chemi- 
schen Struktur der Substanz. 

25. Computerprogrammprodukt nach einem der vorhergehenden Anspruche 19 bis 24, wobei das Vorhersagemodul 
auf einer Datenbasis (108) basiert, die die physikochemischen Eigenschaften von Tests ubstanzen und experi- 

35 mentell fur die Tests ubstanzen ermittelten Parameter beinhaltet. 

26. Computerprogrammprodukt nach einem der vorhergehenden Anspruche 19 bis 25, wobei das Vorhersagemodul 
eine Berechnungsfunktion zur Berechnung des Parameters aus der Lipophilie und/oder dem Molekulargewicht 
der Substanz beinhaltet. 

40 

27. Computerprogrammprodukt nach einem der vorhergehenden Anspruche 19 bis 26, wobei die Berechnungsfunk- 
tion auf einer linearen Regression experimentell ermittelter Parameterwerte basiert. 

45 
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Fig. 2 
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